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We introduce new invariants of a Riemannian singular space, the local 
Yamabe and Sobolev constants, and then go on to prove a general version 
of the Yamabe theorem under that the global Yamabe invariant of the space 
is strictly less than one or the other of these local invariants. This rests on 
a small number of structural assumptions about the space and of the behav- 
ior of the scalar curvature function on its smooth locus. The second half of 
this paper shows how this result applies in the category of smoothly strati- 
fied pseudomanifolds, and we also prove sharp regularity for the solutions 
on these spaces. This sharpens and generalizes the results of Akutagawa and 
Botvinnik Q on the Yamabe problem on spaces with isolated conic singu- 
larities. 

Introduction 

Our aim in this paper is to study a version of the Yamabe problem on a class of com- 
pact Riemannian singular spaces satisfying a small list of general structural axioms 
which we call 'almost smooth metric-measure spaces'. This approach emphasizes 
the centrality of Sobolev inequality, and indeed relies on little else. Our main exis- 
tence result is the analogue of that part of the resolution of this problem on compact 
smooth manifolds (M, g) obtained through the work of Yamabe, Trudinger and 
Aubin, Q, |29l , lf33ll . In that original setting, the work of these authors established 
the existence of a smooth positive function minimizing the Yamabe functional 
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where Scal 9 is the scalar curvature of the metric g, provided the infimum of this 
functional, the so-called Yamabe invariant (sometimes also called the Yamabe con- 
stant or conformal Yamabe invariant) of that conformal class Y(M, [<?]), is strictly 
less than the corresponding invariant of the round sphere. In some papers on this 
subject, the energy Q g is replaced by ((4(n — l)/(n— 2))Q g . The geometric mean- 
ing of this functional is that if u is a minimizer, or indeed any critical point, then 
the conformally related metric g = u n ~ 2 g has constant scalar curvature on any 
open set where u > 0. We refer to the well-known survey paper by Lee and Parker 
ifloll . as well as |[27l . |[28l . for all details on the complete existence theory in the 
setting of smooth compact manifolds. 

The singular spaces (M, g, /x) we are interested in here are typically Rieman- 
nian pseudomanifolds, and in particular Riemannian smoothly stratified spaces 
with iterated edge metrics, endowed with a measure which is a smooth positive 
multiple of the Riemannian volume form dV g . However, as indicated already, we 
require only a few structural assumptions and so our main existence theorem holds 
in much more general settings. In a companion to this paper we explore this di- 
rection further, extending this method to handle rather general semi-Riemannian 
spaces, for example. The ability to allow for a more general measure n is perhaps 
useful, but plays essentially no role in any of the arguments below, and for rea- 
sons of notational simplicity, we often omit [i altogether from the discussion. The 
spaces here are 'mostly smooth' in that they possess an open dense set Q which is 
a smooth ra-dimensional manifold carrying a Riemannian metric. Infinitesimally, 
every point in Q looks the same as every other. However, that is not true if one 
includes the singular points. To accomodate this, we replace the global Yamabe 
invariant by a new invariant which we call the local Yamabe invariant Y^{M, [g]). 
Briefly, this is just the infimum over all points p G M of the Yamabe invariants of 
arbitrarily small balls around p, where we minimize the standard energy functional 
amongst functions on these balls which vanish on the outer boundaries, but not 
necessarily near the singular set of M. We also introduce the corresponding local 
Sobolev invariant Se(M, g). Our main existence theorem states that under various 
sets of conditions on the scalar curvature Scal g (which we regard as a function 
computed in the usual way on the smooth domain if the global Yamabe invari- 
ant Y(M, [g]) is strictly less than Y^(M, [g]) (or, in some versions of the result, 
than Se(M, g)), then Q g admits a strictly positive minimizer u. In certain cases we 
prove that this minimizer u is strictly positive, but show by example that this need 
not be the case if the hypotheses are relaxed. 

Lest this criterion seem too abstract, observe that by conformal invariance, the 
local Yamabe invariant at a smooth point is equal to the Yamabe invariant of the 
round sphere; this is essentially what is known as Aubin's inequality. It is important 
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that Yi involves the limits as r — > of the Yamabe invariants Y(B r (p), g), rather 
than their values at any fixed r > 0; this means that local curvature invariants play 
a smaller role in Yg. An invariant of this nature has been used previously for spaces 
(M, g) with isolated conic singularities. In that setting, if p is a conic point, so that 
some neighbourhood U of p in M is modelled by a cone over a compact smooth 
Riemannian manifold (Z, h), then the local Yamabe invariant at p is the same as 
the so-called cylindrical Yamabe invariant Y(M x Z, [dt 2 + h]) which plays an 
important role in the work of the first author and Botvinnik O, see also (T] for a 
discussion of this problem on orbifolds. It is proved there that Q g has a minimizer 
provided 

- oo < Y(M, [g]) < min{Y(S n , [g ]),Y(R x Z j: [dt 2 + hj])}, (0.2) 
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where (C(Zj), dx 2 + x 2 hj), j = 1, . . . , N, are the local models for the conic 
points of (M,g). We note that Y(R x Zj, [dt 2 + hj}) < Y(S n , [g ]) is always 
true. Also, implicit here is the fact that the cylinder (E x Z, dt 2 + h) and the cone 
(C (Z) , dx 2 + x 2 h) are mutually conformal. We note that there are many examples 
of conic spaces where one does know that d0.2| ) holds, see also IPT51 for some results 
of the case where R g < (rather than Y(M, [g]) < 0), when an existence result is 
obtained in some cases using barriers. 

Our main existence theorem states that if (M, g, fi) is an almost smooth metric- 
measure space which satisfies the first three properties listed in the definition at the 
beginning of §1.1 below as well as one of the three possible hypotheses on Scal 9 , 
then Q g ^ attains its minimum. In certain of these cases, we also prove that the 
minimizing function u is strictly positive on M. The proof is divided into two 
parts: the proof of existence is obtained through a variant of the original method 
appearing in the work of Trudinger and Aubin, and the surprising fact is that this 
original proof may be adapted quite simply to this general setting. However, in 
order to accomodate some of the natural geometric applications later, we present 
an alternate proof of the step which uses Moser iteration to give a uniform upper 
bound for the minimizing sequence, by another argument related to some old ideas 
of Varopoulos. The proof that, in certain cases, the minimizer is strictly positive 
uses some ideas developed by Gursky. 

In the second part of this paper we expand on the theme and setting of Q by 
considering in more detail the case where (M, g) is a smoothly stratified Rieman- 
nian pseudomanifold, also known as an iterated edge space. We identify the local 
Yamabe invariants at all point p £ Mas higher versions of the cylindrical/conic 
Yamabe invariants discussed above; these are simply the global Yamabe invariants 
for the model spaces W k x C(Z), or (conformally) equivalently, M k+l x Z, where 
Z is a compact iterated edge space with lower singular 'depth' than the original 
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space M. The special cases of these invariants when Z = S n ~ k ~ 1 play an inter- 
esting role in the work of Ammann, Dahl and Humbert (5], where quantitative 
estimates of the change of the cr- Yamabe invariant (which is the supremum of the 
Yamabe constants over all conformal classes) under surgeries are obtained. Fi- 
nally, using the more specialized analytic tools available for the study of PDE on 
smoothly stratified spaces, we prove sharp regularity results about the behaviour 
of the minimizer u (or indeed any solution of the Yamabe equation) at the singular 
strata of M. 

The final resolution of the Yamabe problem on smooth manifolds by Schoen, 
described in lfl6l . devolves to showing that Y(M, [g]) < Y(S n , [go]) except when 
(M, g) is conformal to (S n , go). One might hope for some analogue of this result 
here. For example, a natural conjecture is that if (M, g) has only isolated conic sin- 
gularities, then Y(M, [g]) < Yp{M, [g]) unless (M, g) is conformal to the cylinder 
(R x Z,dt 2 + h). Unfortunately, this is now known to be false! Indeed, a re- 
cent paper by Viaclovsky ll32l exhibits a manifold with orbifold singularity which 
does not admit any (incomplete) orbifold metric of constant scalar curvature, or 
equivalently, any finite energy critical point of Q g . Because the existence theory 
in and HI would guarantee a minimizer unless the local and global Yamabe 
invariants are the same, we conclude that these invariants must be equal for the 
spaces Viaclovsky considers. There are no known examples beyond the cylinder 
where Yg(M, [g]) = Y(M, [g]) holds and Q g also has a critical point. It remains a 
tantalizing mystery to determine whether there is some rigidity phenomenon here. 

The authors acknowledge the following grant support: K.A. through the Grant- 
in- Aid for Scientific Research (B), JSPS, No. 24340008; G.C. through the ANR 
grant ACG: ANR-10-BLAN 0105; R.M. through the NSF grant DMS-1105050. 
We are also grateful to Emmanuel Hebey for useful comments. 

1 The general existence theorem 

As discussed above, our main existence theorem yields a minimizer of the func- 
tional Q g on a rather broad class of Riemannian singular spaces. We state and 
prove this result in this section. We first explain the precise geometric and analytic 
hypotheses, then define the local Yamabe and Sobolev invariants and describe their 
relationship to the (global) Sobolev constant of the space in question. We also give 
a number of auxiliary technical facts, including the compactness of the embedding 
W 1 ' 2 L 2p /( p_2 ) for p > n, and that the finiteness of the Sobolev constant im- 
plies discreteness of the spectrum of the (Friedrichs extension of the) Laplacian. 
We also review the standard Moser iteration argument to obtain a uniform upper 
bound for the subcritical solutions and give a different proof based on a different 
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(Morrey-type) assumption on the scalar curvature. Existence of the minimizer then 
follows the lines of the original Trudinger/Aubin argument. The positive lower 
bound for the minimizer uses an argument due to Gursky |[T4l . 

1.1 Almost smooth metric-measure spaces 

Suppose that (M, d, fi) is a compact metric-measure (MM) space which is 'almost 
smooth' in the sense that there is an open dense subset C M which is a smooth 
n-dimensional manifold, and a smooth Riemannian metric g on which induces 
the same metric space structure as d on Q, and hence by density on all of M. (It 
is not hard to check that the arguments in this section only require that g be W 2 ' q 
for some q > n/2, but for simplicity we do not work in this generality.) We also 
assume that the measure dfi is a smooth positive multiple of the volume form, i.e. 
d/j, = h?dV g for some h £ C°°(Q) n C°(M) which is strictly positive. Since g 
and d induce the same distance, we refer to the triple (M, g, fj,) and omit mention 
of d. Note that the metric balls B(p,r) coincide with geodesic balls provided 
B(p,r) C n. 

We shall assume the following properties of the space (M, g, u). 

Hypotheses : 

i) Let W 1,2 (M; dfi) denote the Sobolev space which is the completion of the 
space of Lipschitz function Lip(M) with respect to the usual norm; then we 
assume that Cq($7) is dense in W 1,2 (M; dfi). Notice that this precludes the 
existence of codimension one boundaries. 

ii) Hausdorff n-dimensional measure is absolutely continuous with respect to 
d\i, and both of these measures are Ahlfors n-regular, i.e. 

CrV < fi(B(p,r)) < Cr n 

for some C > 0, and for every p 6 M and r < diam(M). Later, for 
simplicity, we often write fi(B(p, r)) = Yol(p, r). 

iii) The Sobolev inequality holds: there exist A, B > such that 

||/|| 2 _2^<^/ \df\ 2 d^ + B [ \f\ 2 d» (1.1) 

for all / G W 1,2 {M;dn). 

iv) Finally, the scalar curvature function satisfies at least one of the following 
properties: 
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a) Scal 5 6 L 9 (M, d/n) for some g > n/2; 

b) For some q > 1 there exists ana £ [0, 2) such that for every point 

supr" n / | Scal 9 < Cr"" 9 . (1.2) 

r>0 JB(p,r) 

c) Seal" := min{Scal 5 , 0} € L q (M, dp) for some g > n/2; 

We henceforth assume that conditions i) - iii) and at least one of iv) a)-c), are sat- 
isfied, unless explicitly stated otherwise. We call a triplet (M, g, /i) which satisfies 
these properties an almost smooth metric measure space. 

The condition iv) b) states that Scal 9 lies in the Morrey regularity class Ai q x 
where A = n — aq € (0, n]. We state now an important fact about functions which 
lie in this class which will be used in several places below. The proof is deferred 
to the end of §1.3 simply because it involves techniques which are discussed there 
for other reasons, but does not rely on any of the intervening results. 



Lemma 1.1. Suppose that the function V satisfies (1 1 -2b - Then for any e > there 
exists C e > such that for all <j> G W /1 ' 2 (M), 



\V\\ct)\ l dn<e I \d<f)\ 2 dn + C e / \c/)\ 2 dfi. 



Remark 1.2. The extra generality of allowing the measure d[i to be a smooth 
multiple of dV g rather than just the volume form itself plays very little role here. 
For simplicity in this paper, we usually assume that d/j, = dV g . Although the 
analysis in this paper goes through for more general measures, the conclusions 
then are no longer strictly within the realm of conformal geometry. 

1.2 Yamabe and Sobolev constants 

For any open set U C M, we define its Sobolev and Yamabe constants, 

S(U) = inf {f \dip\ 2 dfi : ip <E Wq' 2 (U), \\ip\\ jm_ = 1}, and 

n~ 2 

Y(U) = inf {/(l^p + ^y Seal, ^)dp,: 

if g w 1,2 (wnn),|M|_2n_ = 1}, 

n — 2 

respectively. We also define the local Sobolev constant and local Yamabe invariant 

of (M,g,fj,) by 

S e (M,g) = inf limS(B(p,r)), Y e (M, [g]) = inf Km Y(B(p,r)). 

p£M r->0 p£M r->0 



6 



All these quantities depend on g and u, but we often suppress this, and even explicit 
mention of M, in the notation. 

Lemma 1.3. /fScal 9 satisfies either iv) a) or iv) b), then Y#(M, [g]) = Se(M,g). 

Proof. Assume first that iv) a) holds, i.e. that Scal 5 G L q for some q > n/2. By 
the Holder inequality, 

\Y{U) - S(U)\ < ^^yll Seal, ||, vol(W)^, 

and thus, since q > n/2, 

inf lim S(B(p, r)) = inf lim Y(B(p, r)) 

for any p. 

For the other case, we invoke Lemma [T7T1 as follows. Fix any e > and choose 
C e accordingly. Then 

\Y(B(p,r)-S(B(p,r))\<mi [ Seal, | |0| 2 

<t> JB(p,r) 4 ( n - J-J 

< eS(B(p, r)) + C e inf ( |0| 2 < eS(B(p, r)) + C t 

<t> JB(p,r) 



where the Holder inequality and the normalization ||</>||2n/(n-2) = 1 are used to 
get the last term. Letting r \ and taking the infimum over all p G M shows that 
\Yi — Si\ < eSe, and since this is true for all e > 0, we see that Yf(M) = S%{M), 
as claimed. □ 

The fact that, under these hypotheses, Ye(M) is the same as the local Sobolev 
constant Si(M), leads to an important criterion for the positivity of the local Yam- 
abe invariant. 

Proposition 1.4. Let (M, g) satisfy hypotheses i) - Hi), 
a) For any e > there exists C t > such that 

(Si-e)\\f\\\^ < [ \df\ 2 du + C t f \f\ 2 da (1.3) 
for all f G W 1 ' 2 {M). 
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b) IfSc&lg satisfies either iv) a) or iv) b), then Yi{M) > if and only if the 
Sobolev inequality d 1 - 1 b holds on (M, g, fx). If these conditions are true, 
then for any e > there is a constant C t > such that 

(y £ -e)||/|| 2 _^ < / \df\ 2 dp + C e [ \f\ 2 dp (1.4) 
"- 2 Jq Jn 

for all f £ W 1 ' 2 (M). 

Proof. Let us first address b). If d 1 - lb holds, then rearranging and using the Holder 
inequality, we see that 

S{U) > - (l - B(volU) 2 / n ) , 

which implies a positive lower bound for Yi(M). 

On the other hand, suppose Yf(M) > 0. Fixing any 5 G (0, 1), for each p £ M 
there is a radius r p > such that 

mm{Y(B(p,r p )),S(B(p,r p ))} > (1 - 5)Y e (M). (1.5) 

Since M is compact, there is a finite covering M = [j i B(pi,ri) with rj = ^r Pi . 
Hence if s = min rj, then (11.51) is true for every p £ M with r p replaced by s. 

Now choose a partition of unity {p{\ subordinate to this covering such that 
each Jpi G Lip(M). If / G W 1 ' (M) and x G Lip(M), then x/ G W 1,2 (M). 
Note that this follows from the density assumption i) and the fact that it holds on 
n. Thus for / G W l ' 2 {M), 

(1 - $)Yz(M)\\f\\ 2 2n_ = (1 - ^(M)||/ 2 ||^_ 

< (l-5)Y e (M)^2\\ Pt f 2 \\^_ <^2\\d(^p-f)\\ 2 

i i 

<E(( 1 + Oll^/lli + = (1 + e')|M/||! + C7||/||2 



where (7 depends on the pi and on C e /. This gives (11.41 ) with e = 25Yc(M) if we 
choose e' appropriately. 

The proof of a) is the same. □ 



Remark 1.5. Using the Holder inequality and Lemma [LM we also deduce from 
any one of the hypotheses iv) a), b) or c) that 

Y(M, [g]) > -oo. 
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The Sobolev inequality has another important consequence. 
Proposition 1.6. If the Sobolev inequality holds on (M, g), then the inclusion 

W^ 2 {M) 
is compact for any p G (n, oo). 
Proof. We first write 



u _ e -t(-A+l) u 



(is 



-<~ A+1 Kds 



I e- s (- A+1 \-A + l)uds 
Jo 

f e -| S (-A+l) ( _ A + 1) I e -|s(-A+l) ( _ A + 

Jo 



We shall show that the L 2 ^" 2 ) norm of this difference is bounded by t" \ \u\ | for 
some /3 > 0. This proves that inclusion mapping is approximated in the operator 
norm topology by a sequence of compact mappings, which implies that it must be 
compact. 

To do this, we need three facts. First, if u G W 1,2 (M), then 

(\du\ 2 + M 2 ) = ((-A + l)u,u) = ||(-A + l)^u\\ 2 2 . 
Next, it is an easy consequence of the spectral theorem that 

V-A + le- s (- A+1 ) <Cs-i 

L 2 ->L 2 

Finally, we claim that, if q = 2p/(p — 2) < 2n/(n — 2), then 

C C 



-s(-A+l) 



< 



nt 1_1 ' 



s n/2p ' 



s>0. 



This is proved by interpolation as follows. There is a standard estimate that 



-s(-A+l) 



< 1, 



uniformly in s for any r > 2. This follows by a simpler interpolation from the case 
r = 2 (spectral theorem) and r = oo (easy direct argument). We can obtain the 
same estimate for 1 < r < 2 either by duality or noting that this also holds for 
r = 1 and interpolating again. On the other hand, it is known ll26l that the Sobolev 
inequality (11.41 ) implies that 



-s(-A+l) 



< 



c 

W2' 



s > 0. 
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Thus if we interpolate between this L 1 — > L°° estimate and the U 
with r = 1 — 1/p, then a bit of arithmetic proves the claim. 
Putting these three estimates together, we conclude that 



II estimate 



e -*(- A+1 )u 



< C 



2 2 p ds\\u\ 



1,2 



Ct 2 2 P ||u||l ; 



which decays as required. 



□ 



Corollary 1.7. Let —IS. be the self-adjoint operator obtained as the Friedrichs 
extension from the semi-bounded quadratic form 



(V/,V/) 



M 



over the core domain C^°($7). Then — A has discrete spectrum. 

Proof. It suffices to show that the Friedrichs domain of — A is compactly contained 
in L 2 (M;dp). However, this domain is simply W 2 ' 2 (M; dp) n Wq ,2 (M; dp) C 

W 1,2 (M; dp), which by the previous result is compactly contained in Lp^ 2 for any 
p > n, which in turn continuously includes in L 2 . □ 



1.3 Uniform boundedness of subsolutions 

Let (M, g, dp,) be an almost smooth metric measure space, as considered above. 
We now present two different methods which lead to the uniform boundedness of 
nonnegative functions u which satisfy An > Vu. The first is simply the adaptation 
of the Moser iteration method to this setting, and assumes that V satisfies either 
hypothesis iv) a), or c). This does not cover all the geometric cases we wish to 
consider, so we then prove a stronger result assuming that V satisfies iv) b). This 
second result subsumes the first one, but we describe both proofs since the former 
is the more traditional method and certain constructions in its proof will be used 
later. 



1.3.1 Moser iteration 

We now review the classical Moser iteration method with enough detail to make 
clear that all steps work on almost smooth MM spaces. (In fact, Moser iteration 
works in greater generality still, see [HI.) 

Proposition 1.8. Let u € W 1,2 (M) be nonnegative and satisfy Au — Vu > 0, 
where V G L q for some q > n/2. Then u £ L°° and 

Moo < C(\\V\\ q )\\u\\ 2 , 
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where the constant C depends only on n, q, \ \V\ \ q and the constants A, B from the 
Sobolev inequality. 

Remark 1.9. As usual, the differential inequality is to be interpreted weakly, i.e. 

j (du, dp) <- J Vmp. (1.6) 

for any <p G W 1,2 (M) with p > 0. Notice that the right hand side of (11.61 ) is well 
defined because V G L n l 2 and, from the Sobolev inequality, wp G L n /( n-2 ). 

Proof. We follow the standard proof |fT3l Theorem 8.15] as soon as we verify the 
chain rule: 

Claim: If v G W 1 ' 2 (M) and / G C^R.R) satisfies /' G L°° then / o u G 
fl^M) and 

d(f o u) = f o u . du. 



To prove this claim, note that by [13, Theorem 7.5], we have 

(df o u, dip) = / (du, dip) f o u for all p G Cq(Q), 



and the result follow from the density of Cq(VL) in W 1,2 . 
For a > 2, define 

/.<*>- {*\ - ^ .f»<_f<»-- (L7) 

I x + (a a- 1 — a ) if a Q - x < x; 

the cutoff and additive constant are chosen so that f a (x) is C 1 and convex. Next, 
for any L > 1, let 

a , L (x) = L a f a (f 

Note that (p a ,L(x) = x a on larger and larger intervals as L — >• oo. Furthermore, if 
we define G a ^(x) = J® (j)' a L (t) 2 dt, then a laborious computation gives 

4> a , L (x)<x a <m& xG a , L (x) < -^—-(^(x)) 2 , (1.8) 

2a — 1 

for all x > 0. 
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Inserting ip = G a ^{u) into (11.61) and using (1 1 - 8b . we obtain 

J \d<p a , L {u)\ 2 = J G'^ L {u)\du\ 2 = J(du,d<p) 

< J VuG a , L {u) < ^—j J V<p a , L {uf 



Use both the Sobolev and Holder inequalities to get 

>,2 



,l(u)\\ 2 2„ < 



2 / /■ \ 2=1 



with C = A]| V|L + B(volM) 1 / 9 . Letting L -> oo yields 



n— 2 g-l 

u-j <_^ tt ^j . (1.9) 

This is of course only interesting if the right side is finite. 

We are given the initial choice of q through the potential V, with q > n/2. 
Thus r := 2q/(q — 1) < 2n/(n — 2) and we can choose a sufficiently close to 1 
so that a r < 2n/(n — 2) as well. Since W 1 ' 2 ^ L ar , the right hand side of (11.91 ) 
is finite, and hence so is the left, i.e. u G L Kar where 

n q-l 
K := > 1. 

n — 2 q 



Furthermore, we can rewrite dl.91) as 

IhlUar < (Cia)2^||u|| Q r, (1-10) 

where C\ = (Ca/(2a — 1)). Note that (11.101 ) is valid for any ay > a so long 
as u G L a i T , since aj/(2a>j — 1) is uniformly bounded. Now set ay = K J a, and 
apply (11.101 ) inductively to obtain that 



u 



N-l ^_ 

2kJ 



L"«r < ( II ( C ^ a ) 
j=0 



for any N > 1. Finally, note that the constant here is bounded independently of 
iV; indeed 

tt 1 ,„ ,• \ tt~\ {\ogC\a log^A ^„k — 1 

log n - = £ (-b- + ^rj ^ 

j'=o i=o v 7 
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Thus, taking the limit as N — > oo, we obtain finally that 

IMloo < C{q,n, A, B ,\\V\\ q )\\u\\i2- 

□ 

1.3.2 Varopoulos' method 

For the second method, we now suppose that V satisfies the Morrey condition iv) 
b) for some q > 1 and < a < 2. To compare this hypothesis with the hypothesis 
in Proposition 1 1.81 observe simply that by the Holder inequality, if V £ L n l a , then 
(11.21 ) holds (with Scal 9 replaced by V) provided q < n/a. 

Next, the existence of the Sobolev inequality (11.41 ) implies the Gaussian upper 
bound 

e^(x,y)<C ¥J - 2 e-— (1.11) 

for the Schwartz kernel of the heat operator e tA , see lTT2l . Thus if G(x, y) denotes 
the Green kernel associated to — A + 1, then 

0<G(x,y)< ( 



d n 2 (x,y) 



provided d(x, y) < C. 

Our goal, as before, is to prove the following. 

Theorem 1.10. Assume that u is a nonnegative function in W 1,2 n LP for some 
p > q* such that 

Au > Vu 

where V satisfies (11.21) . Then 

\\u\\oo < C I supsup— / l^l 9 ^) \\u\\ p . 

\ x r >o Vol(x,r) J B(yX ^ J 

(In other words, the constant C depends in some possibly nonlinear way on the 
quantity in parentheses.) 

The key point is to rewrite the differential inequality as (— A + l)u < (— V + 
l)u and then, using that G is positivity preserving, u < G o (— V + l)u. Clearly 
— V + 1 satisfies (11.21 ) if V does, so for simplicity we replace — V + 1 by V. We 
then establish the following mapping properties. 
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Theorem 1.11. Let L, G and V be as above. Then G o V is bounded as a mapping 



n 



L p (M, dp,) -»• L°°(M, dp) when p > 

4 — IX 

L p (M, dp) -»• L»-( 2 - a > (M, d/x), whenpelq*,- , (g* = g/(g-l)) 



2-a 

9 

\ 2 — a 

Proo/ Define the Stieljes measure (i) associated to the nondecreasing function 



i^"x(t)= [ \V(y)\\f(y)\dp(y), 

JB(x,t) 



and then write 

\(G o v)f(x)\ < j M dn ^ y) \ v (y)\ \f(v)\ My) 

where D = diam M. We write the right side of this chain of inequalities as To (/) + 
Too(f), and prove the boundedness properties for these operators separately. 
The estimate of the first of these is trivial. Indeed, for any p > q*, 

MD)<C\\V\\ q \\f\\ p ^-\\T (f)\\ OD <C\\f\\ p . 

We can thus concentrate on T^. 
By the Holder inequality, 

Mr) < Cr~ a \\f\\ p V(x,r) 1 ^ < C\\f\\ p r n ' a -^ . 
Thus if p > 2r^,fhen 

||Too(/)|U = sup [° ^L„ x (r)dr<C'\\f\\ p f ' r 1 ^ dr < C"\\f\\ p . 
x Jo r j 

In the second case, when q* < p < , this integral no longer converges near 
0, so we use instead a classical cutting argument from harmonic analysis, replacing 
the estimate for v x {r) for r small by a different one. 

Define 



v{x)=M(\ffY (x), 
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where M(f) is the maximal function, defined for any L\ oc function / by 

M(f)(x) := SU p / / |/| dp. 
r>0 V(x,r) J B (x,r) 

Hypothesis ii), Ahlfors n-regularity, implies volume doubling, i.e. Vol(x, 2r) < 
C Vol(x, r) for all x € M and r > 0, and from this it is easy to deduce that the 
maximal function, is bounded L l — > £^, cak , 

M{-M(/)>A}<|||/||i. 

Since A4 is also (trivially) bounded L°° — > L°°, by interpolation we see that it is 
also bounded L s — > L s for any 1 < s < oo. This will be invoked below. 
From the definition of v, 

— -a — - 

v x {r) < Cri V(x 1 r)i* v(x) < r n a v(x), 

and hence 



rX(x) Q f D C 

Too(/)(x) < / -^v x {r)dr+ / —^v x {r)dr 

J r J\(x) r 



r\(x) (j pD 
} r Jx(x) 

f\(x) c r D c 

< / — ir —v{x)dr + / 
Jo r Jx 

<C\(xf~ a v{x) + C\(xf- a ~v\\f\\ p 
for any < \(x) < D. The optimal choice of X(x) satisfies 

v(x) = X(x)-^\\f\\ p , 

so inserting this yields 

Too(/)(x) < C\\f\\?- a) " vix) 1 ^ 2 -^ 



n-a-- 



. r n-l' II ^ HP 

AO) ' 



(2-a)r 

pn - 



(1.13) 



p 2 

\Too(f)(x)\*=fc& < C\\f\\^^v{xf. 
Finally, using the boundedness of the maximal function on L s , s = p/q*, we obtain 



\ v \\p ^ ii.' ii/'- 



whence, after some arithmetic, 



ll^oo(/)ll , r , < c \\f\\p- 

n — (2 — a)p 

This is the desired estimate. □ 
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This result implies Theorem 11.101 quite directly. Indeed, reverting back to the 
original potential, we have already noted that < u < Go(- V + l)u and u G LP . 
If p > n/{2 — a), then the first part of Theorem 11.111 bounds | |u| |oo immediately. 
On the other hand, if we only know that p > q*, then the second part of this 
Theorem shows that G o (—V + l)u lies in L Pl where p\ = p (n/(n — (2 — a)p)). 
It is easy to check that there exists e > such that pi/p > 1 + e for any p > q*, 
which means that we can iterate this procedure, obtaining successively that u G L Pj 
for an increasing sequence pj with pj > p(l + e) 3 . Hence p^ > n/(2 — a) for 
some N, so that at the next step u G L°°. 

We conclude this section with the 

Proof of Lemma U~J\ We. begin by noting that under the assumptions of this Lemma, 
the heat kernel bound (11.111 ) holds, hence the Schwartz kernel K fl (x,y) corre- 
sponding to 



In order to prove this Lemma, it is known (see |[25l Theorem X. 18]), that it 
suffices to show that the operator 




for some 5 > 0. Moreover, defining t/? M (r) := Ce jr n 1 , then 




(1.14) 



A„ := (-A + 1 + M 2 ) 1/2 \V\$ : L 2 (X, M ) -> L 2 (X, M ) 



satisfies 



lim ||^4 



= 0. 



Define q* G (1, 2) by ^- + 4r = 1. Then for it G L 2 (X, //), we set 




Because q* G (1, 2), we have that ||<y||2 < C| | tt [| 2- 
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Now introduce the Stieljes measure v x associated to the nondecreasing function 



r4 I |V|2|u|d/i 

JB(x,r) 



so that 

pD i-D 

\(A^u)(x)\ < C / ^(r)cfi? x (r) = C^(D)v x (D) - / </>' (r)^(r) dr, 
Jo Jo 

D = diamM. Using the Holder inequality and the fact that \V\? G L 2 (since 
V G L q for some g > 1), we get 

v x {D) < C\\u\\ 2 , and v x (r) < Cr n ~ a / 2 v{x). 

This gives 

\{A ll u){x)\<C^{D)\\u\\ 2 + [° e-^ r r- a / 2 drv(x) 

Jo 

< C^^D)\\u\\ 2 + C^v{x), 

which proves finally that 

\\AJ L 2^ L2 < CV^D) + <V~f 
Since a < 2, we deduce the result. □ 

1.4 Existence of the minimizer 

We are now in a position to prove the basic existence result. 

Theorem 1.12. Let (M, g, p) be a compact almost smooth metric measure space, 
in particular satisfying hypotheses i) - Hi) in §i.i, and such that Scal 9 satisfies 
either iv) a), b) or c). Supposing that Scal g satisfies iv) a) or b), we then assume 

Y(M, [g])<Y e (M, [g]), 

while if Scalg satisfies iv) c), then our assumption becomes 

Y(M,g)<S e (M,g). 

Then there exists a function u £ W 1,2 {M) n L°°(M) such that ||u||_2n_ = 1 and 



71-2 



Y(M,[g]) = J (\d 



u\ 2 + — ScaL u 2 ) da. 

4(n — 1) 
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Hence on the smooth locus Q, C M, 

n — 2 n — 2 n+2 

An - ^—^ R g u + -^———^Y(M, [g] , ,)u ^ = 



Proof. We follow the lines of the classical proof of Trudinger and Aubin. Since 
W l,2 {M) L 2p ^ p ~ 2 ^ is compact when p > n, the minimum value 

Y v = inf 1 1 (\d<p\ 2 + ^— 2 — Scal 9 ^ dfA 

over all </? £ VF 1,2 (M; d/i) with ||(/?|| 2 2 P = 1 is attained by some function u p . The 



p-2 



usual arguments from the calculus of variations show that u p > and 



n-2 „ . _f±§ 



An p - 77 — ScaL + K, up 3 =0. (1-15) 

4(n — 1) 

It follows from Theorem 11. 101 that u v £ Indeed, dl.151 ) implies that Au p > 
Vu p , where V = c n Scal~ — Y p u 4: ^ p ~ 2 \ Under any of the hypotheses a)-c), Seal" 
satisfies (11.21 ); on the other hand, setting s = p/2q, we have 

r~ n [ \uf~ 2 \ q d^< IKII^r"'/ 5 *-™ =r~ aq 

JB(x,r) p-2 

where a = 2n/p < 2, so the second summand in V satisfies (11.2b as well. Hence 
we may apply this theorem as claimed. 

Now, lim p _j, n l^ = Y{M) < Yi, so for any sufficiently small e > 0, Y p < 
Yi — e provided p is sufficiently close to n. We now may as well replace Y^{M) by 
Si(M), and argue assuming that Scal^ € L 9 for some q > n/2. 

Since ||itp||i,2 is uniformly bounded, we can choose a subsequence — > n 
such that Up. converges to some function u, weakly in W 1,2 and strongly in L q 
for all q G [1, 2n/(n — 2)). Our goal is to show that some further subsequence 



2n 



converges strongly in L n -' 2 . For if this is the case, then we can pass to the limit in 
f n — 2 f 

J {du p , dip) + - Scalg u p cp = Y p J (pu£~ 2 (1.16) 

to conclude that u 6 VF 1,2 n L°° is a weak, and hence strong, solution of the 
equation with ||u||^2_ = 1; setting ip = u into ( 11.16b with p = n then gives that 

n — 2 

Q g (u) = Y(M, [g]), as desired. 
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To accomplish this, recall the function G a ^ introduced in §1.3.1; inserting 
ip = Ga : L(u p ) with a ~ 1 into dl.161 ), and using the properties of these functions, 
gives 



\d(t> a ,L{u p )\ 2 < c(n) J |Scal 9 \\<j) a> L(u p )\ 2 + Y p J u£ +2 \<j) a ,L(u p )\ 2 , 

(1.17) 

where Y p + = max{l^,0} and c(n) = (n — 2)/4(n — 1). Next, using this in the 
Sobolev inequality dl.4| ), with e replaced by e/2, yields that 



(ft - 6/2)11^(^)11^ < / (C e/2 + c(n)| Seal" |) \<P a , L {u P )\< 



a 2 



+ 2^i Y v I < 2 \Mu P )Y 



By the Holder inequality and the normalization of u p , 



p-2 



= \\PaM U p 

n — z 

Furthermore, if \a — 1| and \p — n\ are both sufficiently small, then 

n 2 

— (VolM) 2 ( p -™)/ pn < 1 + e'. 

2a — 1 

Now choose e' so that (1 + e')Y p + < Se — 3e/4. Rearranging the inequality above, 
we obtain 

-rW^a^iup)^^ < C\\4>a yL (u p )\\i + c(n) I I Seal" | \4> a ,L(u P )\ z ■ 

4 n-2 J 

We handle this last term in two different ways, depending on whether Scal 5 
satisfies iv)a) or c), or else iv) b). In the former cases, the Holder inequality es- 
timates this term by || Seal" IWI0a,L(^p)|l2g/( g -iy Since 2a < 2p/(p — 2) and 
aq/(q — 1) < 2p/(p — 2) uniformly as p \ n, we can then pass to a limit as 
L — > oo to conclude that 

< C 

n-2 

for some C which is independent of p. If Scal 9 satisfies iv) b), then we use 
Lemma 11.11 already in (11.171 ) to absorb this term at the expense of an extra fac- 
tor of (1 + e') in front of the Y p + , but we can then proceed exactly as before to 
reach the same conclusion. 
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To conclude, observe that by Proposition (11.81) . and using (11-1 6b . we obtain 
that 1 1 Up | |oo < C with C independent of p. This leads immediately to the strong 
convergence of u p to u in L? n l ( n-2 ) . □ 

Remark 1.13. The existence of a minimizer u E W 1 ' 2 can also be proved using 
the almost optimal Sobolev inequality and a useful trick of Brezis and Lieb ^[9]/. 
However, it still requires the same amount of work as above to prove that this 
minimizer is bounded. 

1.5 The lower bound 

We now show that when Scal 9 satisfies either iv) a) or b), then the minimizer 
obtained in the last subsection is strictly positive. As we show by explicit example 
later, this fails when Scal g only satisfies iv) c). This lower bound is attained by 
adapting an argument of Gursky lfl4l Lemma 4.1]. 

Lemma 1.14. For any ball B(p, r) properly contained in M, there is a constant 
C > such that 

IMU»(ao < c {W\\ L 2 {M) + |M| L2(B) ) \/ip g w l ' 2 (M). 

Proof. Let B' = be the ball with the same center and half the radius. Since 
vol B' > 0, there is a Poincare inequality on the complement of B', 

1Mb < C\\d<p\\ L 2. Vv? E Wq' 2 (M \ B'). 

Hence if p{x) is a Lipschitz cutoff function which equals 1 in B' and vanishes 
outside B, then 

IMIl2(a/) < + IK 1 - p)<p\\l*(m) 

< lklU 2 (B) + C\\d[(l - p)tp]\\i?(M\B>) 

<C (\\d<p\\ L 2 (M) + \\(p\\ L 2 {B) ) 

as claimed. □ 

Proposition 1.15. Assume that Scal 9 satisfies either iv) a) or b). Let u be the 

minimizing solution obtained in the last subsection. Then inf m u > 0. 

Proof. We know that for every <p E W 1,2 (M), 

f f 71 — 2 I n+2 

/ (du,dtp}+ / — — ScaL u<p = Y / u™- 2 (p. 

Jm Jm 4(n - 1) J M 
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It is easy to check that for e, 5 > the function 

cp = ( e + u) 
Inserting this into the identity above gives 



1-28 G W l,2_ 



1 + 25 
and hence 



M 



d(e + u)~ 



n-2 
u 4(n - 1) 



Scal 9 uip + Y I u n ~ 2 ip, 



M 



d(e + u 



-s 



< 



1 + 25 



C 



M 



(| Scal 9 | + u«-2)(e + u) 



M 



-26 



(1.18) 



with C = max j 4 ^_ 2 i) ; 1^1 }• By the Holder inequality, and assuming that Scal s 
satisfies iv) a), 



Ai 



d{e + u)~ 



< 



1 + 25 



a 



26i. 

(e + u) "- 2 



1 1- 



where C" = C [|| Scal 9 || L n/2 + ||«|| i 2n/(«-2)]. If Scal 9 satisfies iv)b) instead, then 
we handle the first term on the right using Lemma [TT1 in an obvious way, and end 
up with the same inequality. 

Applying the Sobolev inequality to the function (u+e)~ s and using the Lemma 
above and this inequality, we conclude that for 5 small enough, 



2Sn 

(e + u) n ~ 2 



<C I (e + u) 

B 



-25 



Assuming that B C O, then by the known upper bound on u and the Harnack 
inequality, there is a c > such that u > c > on 5. Letting e — > in the 
estimate above, we see that 



u »- 2 < C. 



M 



To conclude the proof, note that the convexity of x h- > x 6 implies that 



. n-2 
*4(n- 1) 



ScaJL +5Yu^ u~ d := Vu 



This function V satisfies (11.21 ). so Theorem 11.101 gives that 
inf u > 0. 



< do, i.e. 
□ 
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2 The Yamabe problem on stratified spaces 



We now specialize the results of the last section to the setting of spaces with smooth 
stratifications, also called iterated edge spaces, with corresponding adapted iterated 
edge metrics. We begin by reviewing some aspects of the differential topology and 
metric structure of these spaces, then prove that all such spaces satisfy a Sobolev 
inequality and the other hypotheses from §1.1, and then formulate the precise ex- 
istence theorem in this setting. 

2.1 Smoothly stratified spaces 

We now briefly review the definition of smoothly stratified pseudomanifolds. Fur- 
ther details can be found in in the foundational monograph of Verona PTTl and the 
exposition by Pflaum l24l . Basic definitions vary between sources, and the recent 
paper [4] provides a clarification and unified presentation of some of this material; 
we follow the notation and development of [4, §2] and refer to it for all further 
details, in particular, for a proof that this class of spaces coincides with the class of 
iterated edge spaces considered by Cheeger [ 11 ], cf. also |fl9l . 

Let X be a compact stratified space. By definition, X admits a disjoint de- 
composition into strata, X = USj, where each T>j is a (possibly disconnected, 
possibly open) manifold of dimension j. There are a set of axioms describing how 
the strata fit together, key amongst which is that each connected component of Tjj 
has a tubular neighbourhood U which is the total space of a smooth bundle over 
that component with fibre a truncated cone C(Zj). Here Zj is itself a compact 
stratified space and is called the link of that cone bundle. There is a natural fil- 
tration of X in terms of 'depth' of singularities. Thus compact smooth manifolds 
are said to have depth 0, and if Z is a compact space of depth k, then a space 
which has a neighbourhood which is a truncated cone or a bundle of truncated 
cones with link Z has depth k + 1. This depth filtration is different than the filtra- 
tion of X determined by the closures £j since a stratum of high codimension can 
have low depth (for example, an isolated conic singularity only has depth 1). An 
interesting subtlety is that the fibration of each tubular neighbourhood is required 
to have a smooth trivialization, but it is not a priori obvious what the proper class 
of smooth maps and diffeomorphisms between stratified spaces should be. This 
is precisely the point where the various treatments cited above differ. The defini- 
tion we give is inductive: once a suitable definition of a stratified diffeomorphism 
between spaces of depth j has been given, one declares the suspension of such a 
diffeomorphism, i.e. the radial extension of that diffeomorphism to the cone over 
that space, to be smooth. This extends the definition to spaces of depth j + 1. (This 
smoothness hypothesis excludes many spaces that constitute the standard broader 
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class of stratified spaces, where these local trivializations are only required to be 
continuous; see @ for more on this.) A smoothly stratified space X is a pseudo- 
manifold if the stratum of maximal dimension is dense in X. In distinction to H, 
we allow X to have strata of codimension one, but since hypersurface boundaries 
play a somewhat different role in our main theorem, we say that X is an iterated 
edge space (or smoothly stratified pseudomanifold) with boundary in this case. 

Definition 2.1. For each k > 0, define the class X^ of compact iterated edge spaces 
of depth k as follows: 

• An element of Xq is a compact smooth manifold; 

• A space X lies in X^ if it has a decomposition X = X' U X", where X" is 
an element ofX^_\ with a codimension one boundary along the intersection 
X' n X" and each component ofX' is the total space of a cone bundle over a 
compact base space B with fibre a truncated cone C(Z) for some Z £ Ijt— 1. 
(The common boundary dX' n dX" is the total space of a bundle over the 
same base B with fibre Z.) 

• Any element iGlt has a well-defined dimension, where in the decomposi- 
tion above, dim X" = dim X' = dim B + dim Z + 1. 

Note that if X £ then its stratum Y of maximal depth k is necessarily a 
compact smooth manifold. 

Every iterated edge space X carries a class of adapted iterated edge metrics, 
which are also defined inductively. Thus assuming that we have described the 
class of admissible iterated edge metrics on all iterated edge spaces of depth k — 1, 
let X be an iterated edge space of depth k. If Y is the stratum of depth k and 
U the tubular neighbourhood around Y, then we can assume that the structure of 
the metric g on X has been described on X \ U. In particular, if x is the radial 
coordinate on the conic fibres of U, then dU = {x = 1} is an iterated edge space 
of depth k — 1, which is the total space of a fibration over Y with fibre Z. Let G be 
an admissible metric on dU, which we assume has been defined by the inductive 
hypothesis. Then it is of the form w*h + k, where it : dU — > Y is the fibration, h is 
an ordinary Riemannian metric on Y and k is a symmetric 2-tensor on dU which 
restricts to an admissible metric on each fibre Z. We also assume that k is totally 
degenerate on a subspace of dimension I at each point, I = dim Y. Now define the 
metric g on U by coning off each fibre. In other words, we set 

<7o = dx 2 + n*h + x 2 k, 

where x is the radial function on each conical fibre. An admissible metric on U 
is any metric which has the form g$ + n where k is polyhomogeneous on the 
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resolution of this neighbourhood and such that \k\ 90 decays at some rate x 1 . The 
notion of polyhomogeneous regularity will be defined in the next section, so for the 
moment consider this only as an appropriate smoothness condition. It is possible to 
consider finite regularity metrics too, but for simplicity we shall not do so. Finally, 
a metric g on all of X is admissible if it is admissible in a sense defined via the 
inductive hypothesis away from hi, and which takes this form in hi. 

To make this more explicit, let V x Z be a local trivialization of dU, where 
V C M. e is an open ball, and Z is the depth k — 1 link, and introduce a coordinate 
system y G V as well as local coordinates z on the smooth stratum of Z. Then we 
can write 

l l n-e-l n-e-1 

G = Yl h ij(y)dy i dy j + Y, £ b ip (y,z)dy i dz p + ^ k pq (y, z)dz p dz q , 

i,j=l i=l p=l p,q=l 

where the second and third sums here constitute the tensor k. Thus 

I 

gi = dx 2 + ^2 h ij (x,y)dy l dy J + 

i,j=l 

l n-e-i n-e-i 
x<2 /^2 t } ip(x,y,z)dy t dz p + x 2 ^ k pq (x,y, z)dz p dz q , 
i=i p=i P)<?=i 

where hij, bi v and k pq are smooth for < x < 1. 

Pick any point p G Y, the depth k stratum, and use coordinates so that x = 0, 
y = and z = zq at p. Define the dilations D\ : (x, y, z) — > (Xx, Xy, z), and 
consider the family of metrics g\ := Dy x g\. Then, as A — > oo, 

I n-e-i 
A 2 5a — »• ^ 2 + Yl h ij {0)dy i dy j +x 2 k pq {0, z)dz p dz q 

i,j=l P,q=l 

which is a product metric on the space R e x C(Z) (where M. e is identified with 
T p Y). Note that the metric k on Z in this product decomposition depends on the 
basepoint p G Y, whereas h(p) is simply the Euclidean metric in some linear 
change of coordinates. We summarize this by saying that \dy\ 2 + k(y) is the model 
iterated edge metric for g at y G Y. Observe that the perturbation k disappears 
in this same rescaling limit. An important consequence of this metric structure is 
that we can choose local coordinates (x, y, z) near any fibre Z y corresponding to 
y G Y such that the scalar Laplacian takes the form 

A g = dl + x" 1 A(x, y, z)d x + x' 2 A k(Xty ^ z) + A h (,x, y ,z) + E , 
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where all coefficients are smooth, or at least bounded and polyhomogeneous, with 
A(0, y,z) =n — l—l, and where E is a higher order error term of first order in 
the sense that it is a sum of smooth multiples of the vector fields xd x , xd y and d z . 

There are slightly less restrictive types of metrics which one can handle without 
too much more difficulty; for example, one could allow terms like xdy l dz p , or (for 
the final metric, after the perturbation k is added), terms like dxdy 1 or xdxdz p , but 
again for simplicity we do not do so here. 

2.2 Sobolev inequalities 

We next show that the Sobolev inequality (11.11 ) holds on any iterated edge space 
with adapted metric. 

Proposition 2.2. Let (M, g) be an iterated edge space, possibly with boundary, 
with admissible metric g as defined in the last subsection. Denote by Q, its principal 
open dense stratum. Then the Sobolev inequality (11.11 ) is valid for all u 6 Qf°(f2), 
and hence for all u 6 Wq' 2 (M). 

Proof. We reduce the problem of verifying (11.11) on an iterated edge space (M, g) 
of depth k using the following observations. First, the Sobolev inequality is local- 
izable; in other words, if (11.11) holds on every set in a finite open cover {U a } of M, 
then using a partition of unity we can show that it holds on all of M. Now decom- 
pose M = M'uU where U is the tubular neighbourhood around the maximal depth 
stratum in M and M' is an iterated edge space of depth k — 1 with boundary. We 
may assume by induction that (11.11 ) holds for all functions with support in M', so it 
suffices to verify this inequality for functions with support inU. Localizing further, 
we can restrict attention to functions supported in a local trivialization V x C\(Z) 
of U, where V C M. e is an open ball and Z is a compact space of depth strictly less 
than k. Finally, noting that (11.11) is stable under quasi-isometric changes of metric, 
we may assume that g is the product metric \dy\ 2 + dx 2 + x 2 hz on M. e x C(Z). 

We now recall the fact that (11.11 ) holds on a space (W, gw) if and only if the 
heat kernel H w (t, w,w') for the scalar Laplacian satisfies 

H w (t,w,w') <C'r n ' 2 (2.1) 

for all w,w' G W and < t < 1, where n = dimW. (Indeed, [26, Theorem 
4.1.3] states that the Nash inequality is equivalent to this heat kernel estimate; the 
equivalence of the Nash inequality with (11.11 ) is treated in |26l Ch. 3] ; alternately, 
ll2ll shows that the Sobolev inequality implies the heat kernel bound, while by 
|[30ll . the heat kernel bound implies the Sobolev inequality.) 
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We apply this in two separate ways. First, since Z is a compact iterated edge 
space of depth less than k, dl.ll ) holds on Z; hence H z (t, z, z') < Ct~ m / 2 where 
m = dim Z = n — I — 1. Using this, we shall show that 

H c( ~ z \t,x,z,x',z') < Ct^ m+1 ^ 2 . (2.2) 

Since the corresponding heat kernel bound on W is standard, and since heat kernels 
multiply for Riemannian products, we see that 

H^ c ^ z \t,y,x,z,y',x',z') = 

H Re (t,y,y')H c ^(t,x,z,x',z') < Cf'/Mn-fl/a = Ct~ n / 2 . 

Hence O) holds on V x Ci(Z). 

It remains to verify (12.2b - Denote by H a,b the heat kernel on the conic nappe 
C a fi(Z) = {(x, z) : a < x < b}, with Dirichlet conditions at the boundaries. Note 
that (11.11 ) holds on C\^{Z) with respect to the product metric, hence by quasi- 
isometry invariance, it also holds with respect to the conic metric. Therefore, 

H 1 ' 2 ^ x, z, x', z') < Ct^ m+1 ^ 2 . (2.3) 

Now recall the basic scaling property of the heat kernel. For any A > 0, the heat 
kernels H X,2X and H 1,2 are related to one another by 

^ A ' 2A (A 2 t, Ax, y, Ax', y')\ m+l = H l,2 (t, x, y, x' , y'). 

Using (12.31 ) and changing variables, we obtain 

H x > 2X (t,x,y,x',y') <Ct^ m+1 ^ 2 . 

The squared L 2n /( n ~ 2 ^ norm on the left in (11.11 ) and the squared L 2 norm of u on 
the right both scale the same way, but the squared L 2 norm of Vn scales differently. 
Thus when we apply this for the sequence A = 2 _J , and assemble the pieces using 
a dyadic partition of unity {x(2- J x)}, where x i s supported on 1/2 < x < 4, then 
we conclude that 

n-2 



(/ x- 2 u^dV g ) <C\[ \Vu\ 2 dV g + [ x- 2 u 2 dV g ), 

\JC ,i(Z) J \JCo,i(Z) JC 0rl (Z) J 

which is valid for all u £ C^(C 0A (Z) \ {0}). 

fn— 2) In 

Since x < 1 in the support, the left side dominates 11^1 l2n/(ri— 2) • ^ n me omer 
hand, we claim that there is a Poincare-Hardy inequality in this setting, i.e. 



(m- if 



x- 2 u 2 dV g < I \Vu\ 2 dV g . 
C(Z) Jc(z) 
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This is standard when Z is a compact smooth manifold, but since Az has discrete 
spectrum by virtue of Corollary 11.71 and the inductive hypothesis, we can reduce 
to the individual eigenspaces, where it becomes the usual Hardy inequality on M + . 
As an alternate path to proving this we could use the argument in |[T0l , which uses 
integration by parts and hence requires only the density of functions with compact 
support in the smooth locus. 

We have now verified (12.21 ). and hence have proved that (11.11 ) holds for all 
iterated edge spaces of depth k. □ 

Appealing to Lemma [2~4l below and combining the result above with Proposi- 
tion [L4j we obtain the 

Corollary 2.3. Let (M, g) be a compact iterated edge space with adapted metric 
which satisfies one of the conditions in Lemma \2.4\ so that at least one of the hy- 
potheses iv) a) or iv)b) hold. Then the local Yamabe constant Y^{M, [g]) is strictly 
positive. 

In the next subsection we identify this local Yamabe constant somewhat more 
explicitly. 

2.3 Existence of Yamabe metrics 

We now turn to the problem of finding minimizers for the functional Q g in this 
setting of iterated edge spaces. The main issue now is to understand when the 
hypotheses iv) a), b) or c) hold so that we can apply Theorem 1 1.121 

We first describe the local Yamabe invariant of an iterated edge space (M, g). 
Let p G M. If p lies in the depth stratum, i.e. is a smooth point, then the local 
Yamabe invariant at p is just Y(S n ). If p lies on a depth k stratum E, then as 
described at the end of §2.1, the rescaled limit of the metric g equals dx 2 + dy 2 + 
x 2 k p , where k v is the metric on the link Z at p and dy 2 is the Euclidean metric 
on IR , t = dimE. Note that this is conformally equivalent to the product metric 
g m t+i + k p on H^ +1 x Z, and hence 

Y(R e x C{Z), dx 2 + dy 2 + x 2 k p ) = Y(H m x Z, [g m t+i + kp]). 

This generalizes the fact that the Yamabe invariant of the cone C(Z) and the cylin- 
der M x Z are the same. In any case, enumerating the depth j strata as and 
denoting the link around Sj by Zj, then we have proved that 

Y e (M,g) = min inf {Y(R e x C(Zj), [dy 2 + dx 2 + x 2 (k z X])}. (2.4) 

Now consider the hypotheses in § 1 . 1 . The verification of i) is a straightforward 
exercise using cutoff functions and mollifications, which we leave to the reader. 
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The Ahlfors ra-regularity is even easier. We have verified in §2.2 that the Sobolev 
inequality dl.ll ) holds; this is condition iii). On the other hand, the hypotheses iv) 
a)-c) require more careful attention. Indeed, as we now show, these hypotheses are 
valid for a rather limited set of iterated edge metrics. 

We begin with some general remarks. It is clear from the structure of adapted 
iterated edge metrics that if Yij is any stratum of depth j and xj is the radial distance 
function in the tubular neighbourhood of Ej, then the scalar curvature Scal 9 can 
blow up no faster than xj 2 . If dim Zj = fj, so dim Sj = tj = n — fj — 1, then 

dV g « Xj dxjdVh-dV^j near Sj, where hj is a smooth metric on T.j pulled back 
to the tubular neighbourhood and kj restricts to a metric on the (depth j — 1) link 
Zj . Assuming that g is smooth in the variable Xj , then 

A (J) A U) 
Scal g = ^ + ^ + 0(l). 

Xj Xj 

To correlate this with the hypotheses iv) a) - c), note that € L q implies 
q < (fj + l)/2, and hence we can never take q > n/2 as in iv) a). Similarly, the 
Morrey condition requires that for some q > 1, 

r~ n j xj 2q+f] dxjdydz = Cr - n+f i- 2q+l+l i = Cr~ 2q , 

which is (11.21) with a = 2 and hence does not fit into our hypotheses. Suppose, 
however, that the coefficient of xj 2 in this expansion vanishes. Then xj 1 € L q 
provided q < fj + 1, and hence we can take q > n/2 and see that iv) a) is satisfied 
provided fj + 1 > n/2. Similarly, the Morrey condition holds because we only 
need choose q > 1, which is always possible since fj + 1 > 1, and for such a q we 
then have 

r-' n f xj q+f] dxjdydz = Cr- n+f i~ q+l+ ^ = Cr- q , 
Jb,. 

which is (11.21 ) with a = 1. This proves the 

Lemma 2.4. The scalar curvature Scal 9 satisfies iv) a) if and only if then = 
for all j and in addition A^p = whenever fj < (n — 2)/2. The scalar curvature 
Scal 9 satisfies iv) b) (for some q > 1 and < a < 2) if and only if = Ofor 
all j. Finally, Scal 9 satisfies iv) c) if and only if Aq^ > Ofor all j and A^ > 
when fj < (n - 2) /2. 

It is clear from this that the terms in Scal 9 which blow up like are the 
most problematic. The l/xj terms always fit within hypothesis iv) b). 
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We will also need to consider metrics with a polyhomogeneous expansion, 
which include noninteger powers of x or terms like x 7 (logx)^, I 6 No- For any 
such metric, the scalar curvature function also has an expansion and there is an 
obvious extension of this lemma which requires the vanishing or nonnegativity of 
the coefficient of any term x 7 (log x) e where 717/2 + fj < 1. 

Isolated conic points 

We first examine the simplest case: an isolated conic singularity, where dim Sj = 
and fj = n. For simplicity drop the index j, but to be consistent with later 
notation, we still use / = n — 1. A well-known formula 11231 p. 69] shows that an 
exact warped product conic metric g = dx 2 + x 2 k has Scal 9 = x~ 2 (Scalfc — /(/ — 
1)). More generally, if k depends smoothly on x, then 

Scakfm — f(f — I) , 

Seal, = () , + Oix- 1 ) (2.5) 

x l 

This leading coefficient vanishes if and only if Scal/^o) = /(/ — 1), which indi- 
cates a very strong geometric and topological obstruction: if the scalar curvature of 
(M, g) is bounded, then in particular the link (Z, k(0)) must have positive Yamabe 
invariant. 

We now study whether it is possible to remove the singular terms in the expan- 

^ 4 

sion of Scal 9 using a conformal change. If g = w n ~ 2 g, then 

1 n + 2 77, — 2 

Seal, = — c(n) w ™- 2 (A„w — c(n) ScaL w), c(n) = — -. (2.6) 

4(n — 1) 

Thus if we introduce the expansion in x of A g and Scal ff , we obtain that 

Scalg ~ — c(n)~ w n ~ 2 x 

(d 2 x + —d x + ^(A ko - c(n)(Scal fc(0 ) -/(/ - 1))) + -e) w. 

V X X X J 

The error term E ~ Eq+xE\-\- . . . is a second order differential operator composed 
of a sum of smooth multiples of products of the vector fields xd x and d z , and also 
includes the terms beyond the leading one in the expansion for Scal ff . 

From this we see that a necessary and sufficient condition for the coefficient of 
x~ 2 to vanish is that 

(Afc(o) - c(n) Scal fe( o))u;o = -c(n)/(/ - l)w , 
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where wq is the restriction of w to x = 0, i.e. the leading term in the expansion of 
w. We denote the operator which appears on the left here by £%^ Q y it is a special 
element of the family of operators 

Tfl — 2 

£™o = A ko - c(m) Scal fc(0) , c(m) = ^ _ -r , (2.7) 

for any value of to. Note that £L, is simply the conformal Laplacian of (Z, k(0)). 
The positivity of the operator — j£j^ ) plays an important role in the main existence 
theorem of 0, as we now recall. 

Theorem 2.5 ([3 ]). Suppose that (M n , g) is a space with isolated conic singular- 
ities, and that at each conic point p, the operator — ^( ) on ^ e M^)) 
a// eigenvalues strictly positive. Suppose too that Y(M, [g]) < Yi(M, [g]). Then 
there exists a function uon M which minimizes Q g and is such that g = u n ~ 2 g re- 
mains incomplete. Conversely, there exists a minimizer u such that g is incomplete 
only if-£ k (o) > 0. 

We do not assert that u is bounded, nor that the new constant scalar curvature 
metric is conic. We shall explain shortly why u may fail to be bounded; in the next 
section we describe the polyhomogeneous regularity of u which makes clear that g 
is in fact still conic. 

We can recover part of this theorem immediately from Theorem 11.121 Indeed, 
if the lowest eigenvalue of —^(o) ^ s exact h c ( n )f{f ~ 1). then we can choose 
w so that wq is the eigenfunction corresponding to this lowest eigenvalue, so that 
wq is strictly positive, and then Scal^ blows up only like x^ 1 , hence lies in L q for 
q G (n/2, n). We thus obtain the existence of a bounded, strictly positive function 
u which minimizes Qg. Clearly u*'( n ~ 2 'g is quasi-isometric to g and thence to g. 

In order to prove existence whenever — £w Q x is positive, fix 5 > 0, to be spec- 
ified below, and define g§ = x 2S ~ 2 g. The change of variables £ = x s /5 gives the 
transform 

gs = x 25 - 2 (dx 2 + x 2 k) = df + f5 2 k, 

so gs is still conic, but its link metric has been scaled by 5 2 . Oserve also that 
— Cfi2 k t \ = — ^ _2 ^fc(o)- This means that if we first replace g by gs and then set 
gs = u; 4 /( n-2 ) gs, then we can make the coefficient of vanish provided that 
5~~ 2 \o(— jCwqj) = c(n)f(f — 1), which determines the value of 5, and wo is the 
corresponding eigenfunction. We are then in a position to apply Theorem 11.121 
again, this time with gs as the background metric. The solution metric is quasi- 
isometric to gs, and hence conic. 
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The converse statement is an easy consequence of these same calculations, at 
least once we show that a minimizer u (or indeed any positive solution of the corre- 
sponding Euler-Lagrange equation) has a polyhomogeneous expansion as x — > 0, 
which we do in the next section. Thus we have now given an independent proof of 
Theorem 12.5 1 

The condition that Ao(— ££/ p > is actually stronger than the condition that 
(Z, k(0)) is Yamabe positive. Indeed, referring back to the family of operators 
(12.7I ). an easy calculation shows that if p < q, then there are positive constants 
A = A(p, q) and B = B(p, q) such that 

" ^(o) = ^(-4(0)) + B(-A k{0) ) =► -£* m > A{-C\ m ). (2.8) 

In particular, taking p = f and q = n, then the positivity of — £^( ) implies that 
the conformal Laplacian of (Z, k(0)) is positive, which is well-known to imply the 
existence of a conformally equivalent (constant) positive scalar curvature metric. 

Simple edges 

We next suppose that M has only simple edges, i.e. that each singular stratum Sj 
is a compact smooth manifold of dimension n — rj. For simplicity we assume 
that there is only one such stratum and drop the index j. A tubular neighbourhood 
of £ is a cone bundle with compact smooth link Z* , and in this neighbourhood, 
g ~ dx 2 + x 2 k + -K*h, where h is a metric on X, it the projection from this 
neighbourhood onto S and k a symmetric 2-tensor so that dx 2 + x 2 k pulls back to 
an asymptotically conic metric on each conical fibre of the tubular neighbourhood. 

Lemma 2.6. If g has a smooth expansion as x — >• 0, then 

Seal, = SCal ^°^ : /(/ " 1} + + O(l). (2.9) 

x z x 



This is slightly less obvious than in the isolated conic case and can be veri- 
fied by direct calculation. It can also be proved by observing that since Scal 9 has 
an expansion with initial term x~ 2 , if we dilate the coordinates via (x, y, z) — > 
(Ax, Xy, z) (around some fixed basepoint y$ G S) and let g\ be the corresponding 
pulled back metric, then the coefficient of the leading term of homogeneity —2 in 
A must be the limit of A 2 Scal 9A as A — > oo. However, it is evident that \~ 2 g\ 
converges to the product metric dx 2 + x 2 k(y , z) + dy 2 on C(Z) x M n_r , which 
has scalar curvature exactly equal to x~ 2 (Scalfc( ) — /(/ — 1))- 

Let us now investigate whether it is possible to conformally transform away the 
singular term of order 1/x 2 in the expansion for Scal 9 at S. As we have already 
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shown, the existence of the singular term A\jx can be handled using hypothesis 
iv)b). 

^ 4 

Replace g by g = w n ~ 2 g and proceed with exactly the same formal calcula- 
tion as in the isolated conic case. There are several important differences in this 
setting. First, it is still clearly necessary that Ao(— £&(o)( y )) = c ( n )f(f ~ 1) an d 
that wo(y, z) must lie in this eigenspace for every y. In particular, this eigenvalue 
must be independent of y £ S, which is a strong rigidity statement. Assuming this, 
we can thus eliminate the x~ 2 term. Note that we may - and indeed we shall later 
need to - let wq depend nontrivially but smoothly on y. This does not interfere with 
this calculation since although E now contains y derivatives, these are accompa- 
nied by a nonnegative power of x, hence the derivatives of wq can be regarded as 
junk terms in the expansion and can be solved away. 

Applying Theorem [TTT2] and the fact that we have arranged that Scal 9 satisfies 
iv) b), we have now proved the 

Theorem 2.7. Let (M,g) have at most simple edge singularities. Assume that 
c(n)f( f — 1) = Ao(— £fc .(o)) along each singular stratum Sj. Suppose in addition 
that Y(M, [g]) < Yi(M, [g]). Then there exists a bounded and strictly positive 
function u which minimizes Qg. The metric ii 4 '( n ~ 2 )<? is quasi-isometric to the 
initial metric g. 

Unlike the conic case, we cannot go further and still remain within the class 
of iterated edge metrics. Indeed, if we were to multiply g by the conformal factor 
x 25 ~ 2 , then this factor would also multiply ir*h; if 5 < 1, the corresponding metric 
would have infinite diameter, while if 5 > 1 then the entire edge £ would be 
collapsed to a point. In either case, we would leave the category of smoothly 
stratified spaces and iterated edge metrics. 

We shall not carry out the detailed study of when we can modify g conformally 
to ensure the weaker condition iv) c), that (Scal^)_ G L q for some q > n/2. The 
conditions are not particularly explicit, and the solution u is not bounded away 
from so that the solution metric is again not of iterated edge type. 

The general case 

We now come to the general case where (M, g) is a smoothly stratified space with 
iterated edge metric. As we shall explain, the conditions needed to obtain a solution 
of the Yamabe problem in this category are even more restrictive than in the simple 
edge case. 

We begin with a statement of the simplest case, which follows immediately 
from Theorem 1 1.121 
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Theorem 2.8. Let (M,g) be a compact smoothly stratified space with iterated 
edge metric g. Suppose that along each stratum £,-, the link metric (Zj, kj) has 
Scalfe- = fj(fj — 1), and in addition, that Y(M, [g]) < Ye(M, [g]). Then there 

4 

exists a bounded, strictly positive function u which minimizes Q g , and hence u n ~ 2 g 
is an iterated edge metric with constant scalar curvature. 

The regularity theorem in the next section will show that u is polyhomoge- 
neous, so that this solution metric is indeed an iterated edge metric in the strict 
sense of the word. 

As in the conic and simple edge cases, we might also seek conditions on the 
initial metric g so that there is some conformally related metric g which satisfies the 
hypotheses of Theorem ll.121 As in those cases, the idea is to choose the conformal 
factor w to kill the appropriate singular terms at each stratum. 

The calculations we have done above may be carried out almost exactly as 
before, and lead to the following necessary conditions: for any stratum £ with link 
(Z? , k), we assume that 

1) The operator — £^ on Z has discrete spectrum; 

2) The operator — £J? has lowest eigenvalue c(n)/(/ — 1) at every point of S; 

The first hypothesis, on the discreteness of the spectrum, may be surprising. The 
fact that the scalar Laplacian itself is essentially self-adjoint and has discrete 
spectrum is a consequence of Corollary 11.71 and Proposition 12.21 However, the 
extra term c(n) Scalfc may blow up like 1/r 2 on approach to any of the singular 
strata of Z itself, which changes the indicial roots. It is not hard to find examples 
of spaces (Z, k), even with just isolated conic singularities, where —C^ is not even 
semi-bounded, which simply amounts to the fact that c(n) Scal^ diverges to — oo 
like —c/r 2 with leading coefficient larger c than the permissible Hardy estimate 
bound (/ — l) 2 /4. This question is closely related to the problems studied in [20], 
see also IM and 

One further point which requires explanation is that in using condition 2), we 
use a conformal factor w which has leading coefficient along £ equal to the eigen- 
function wq for — corresponding to the eigenvalue c(n)f(f — 1). In order to stay 
with the class of iterated edge metrics, it is necessary that wo be bounded above 
and strictly positive, and this may fail. Indeed, it is easy to construct examples of 
operators — + V on Z with V blowing up like 1/r 2 , where the ground state 
eigenfunction either vanishes at the singular set of Z or else blows up at some rate. 
Fortunately, the fact that this does not occur follows from the hypotheses we have 
already made. 
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Proposition 2.9. If (M,g) satisfies conditions I) - 2) along each singular stra- 
tum E, and if (Z, k) is any link, then the eigenfunction wq for the ground state 
eigenvalue of —C^ is bounded and strictly positive. 

Proof. The assertion follows from the regularity theory for eigenf unctions, re- 
viewed in the next section, and an indicial root computation. Let £' be any sin- 
gular stratum of Z with corresponding link (Z', k'), dim Z' = f. If r is the radial 
variable to this stratum, then near S', 

C n k = d 2 r + ^-d r + l(Afc, - c(n)(ScaW -/'(/' - 1))) + A E , + E' , 

where E' contains all higher order terms (including higher order terms in the ex- 
pansion of Scalfe). The indicial roots of this operator are then equal to 

± /'-i / (/'-i) 2 , 

v i = — — ± y^^ + ^ 

where the fij are the eigenvalues of — (A^/ — c(n)(Scalfc/ —/'(/' — 1)) = — — 
c(n)f'(f — 1). By assumption, /Uo = < < • • • , hence 

i/jj" = < < . . . , and = 1 - f > z/{~ > — 

By the aforementioned regularity theory, wo ~ cr u ^ + (j)j + . . . near S', where 0j is 
the eigenfunction corresponding to i/T. However, must remain strictly positive 
in the interior of Z by the standard maximum principle arguments, hence j = 
and wo ~ c^o which shows that it remains bounded and strictly positive near this 
stratum. □ 

To conclude this section we observe finally that assuming the conditions 1) and 
2) on (M , g), if (Z, k) is any link, then by (12.8I ). the conformal Laplacian —C k is 
strictly positive. 



3 Regularity 

The final goal of this paper is to study the regularity of the minimizers of the 
functional Q g obtained in the last section when (M, g) is an iterated edge space. 
The techniques here are nonvariational, so the results below apply to any positive 
solution of 

n + 2 

A g u — c(n) Scal 9 u + c(n)Au™~ 2 = 0, (3.1) 

4 

assuming that u satisfies a natural growth condition so that u n ~ 2 g remains quasi- 
isometric to g, and which is satisfied for the solutions constructed in §2.3. Note that 
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if ||w||_2n_ = 1, then A = Scal 9 > Y(M, [go]). We shall prove that u is conormal 
along each of the singular strata, and has (at least) a partially polyhomogeneous 
expansion. We explain this below. 

As in §2.3, we first prove regularity when M has only isolated conic singular- 
ities. The steps in this case are quite elementary, but rely on a certain number of 
definitions concerning the function spaces and the 6-calculus of pseudodifferential 
operators. With these preliminaries, the proof of regularity in this case is only a 
few lines. We then prove regularity for spaces with simple edges, and here we can 
quote known results about the pseudodifferential edge calculus from iMTl . For the 
general case we need only mimic one small part of this edge calculus to be able to 
deduce what we need. 

3.1 Conic singularities 

Our first goal is to prove the 

Proposition 3.1. Let (M, g) be a compact space with only isolated conic singu- 
larities. Assume that g is a polyhomogeneous conic metric. Suppose that u is 
a solution of (13.11) which is positive on the regular part of M and which satis- 
fies u < Cx~( n ~ 2 >' 2+€ for some e > near each conic point, where x is the 
radial distance to the conic tip. Then u is polyhomogeneous as x —> 0. If the 
link (Z,k) satisfies the simplest condition, that = (n — l)(n — 2), then the 
expansion of u takes the form u ~ cq + c\(z)x Vl + . . ., where cq is a positive 
constant. If we only have that the lowest eigenvalue of —£ r £ is positive, then 
u ~ x^~^ n ~ 2 ^' 2 [co{z) + ci(z)x u 't +...), where cq(z) is strictly positive and 
is the ground state eigenfunction of —C^ and 6 is the constant described in §2.3. 
The exponents ui, u[, etc., which appear in this expansion are determined by the 
higher eigenvalues of 

This regularity is local near each conic tip, but we emphasize that it is global 
with respect to the links (Z,k). The main issue is to prove that the solution is 
conormal (see below); its polyhomogeneity and the precise form of its expansion 
are then formal consequences. 

Rather than analyzing (13.11 ) directly, we rewrite it relative to the background 
metric g = x~ 2 g = (a^ 2- ™)/ 2 ) 4 /*™ -2 ), yielding 

(Ag - c(n)Rg)v + c(n)Aw«-2 = 0, (3.2) 

where the original solution u = x^ 2 ~ n ^ 2 v. Note that by the transformation prop- 
erties of the conformal Laplacian, Rg = Ruq\ + • • • , hence 

((x^) 2 + A fc(0) - c{n)R m +E)v + c{n)Av^ = 0; (3.3) 
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here E contains all higher order error terms from both and Scal^. Note that 
our assumption that \u\ < C'x^ 2 ~ n ^ 2+t becomes \v\ < Cx e , which is much easier 
to work with. (In fact, it is precisely because we are able to work with solutions 
which decay that this argument is easier than the corresponding regularity theorem 
in HI.) 

Before embarking on the proof, we recall several facts, first about the function 
spaces which will be used and then about parametrices in the 6-calculus. For sim- 
plicity, assume that M has only one conic point and that the radial function x is 
extended globally and is strictly positive elsewhere on M. 

Definition 3.2. Decompose M as M' U C\(Z), where the second factor is the 
truncated cone over Z with coordinates z G Z and x G (0, 1]. It is most natural to 
work relative to the complete metric g, and in this geometry, — log x is the distance 
function on the asymptotically cylindrical end. 

i) The space C^' 7 (M) consists of all functions v which lie in the ordinary 
Holder space C fc ' 7 on M', and in addition satisfy (xd x y d^v G C^' 7 , where 
the latter space is defined using the seminorm 

\v(x, z) — v(x', z')\ 
(x,z)^(x>,z>) distg((x, z), (x 1 , z')) 7 ' 

l/2<x/x'<2 

We also define 

x^(M) = {v = x^v : v E C£' 7 (M)}. 

ii) For any u E R, let A U (M) = f] k>0 x u C k '^(M). This is the space of 
conormal functions. Next, define the space of polyhomogeneous functions 
^4phg(M) to consist of all conormal functions v which admit complete asymp- 
totic expansions with smooth coefficients, and write Ap hg for all polyhomo- 
geneous functions with leading term x u °4>(z) for some vq with (real part) 
greater than or equal to v. Note that x v log x G -^phg f or an y e > 0. 
Finally, let v < v' he any pair of real numbers, and define -4p^g C^O = 
.4p h (Af) + A u ' (M); thus v is in this space if it has a partial polyhomo- 
geneous expansion with initial term bounded by x v and with conormal 're- 
mainder' vanishing like x u . 

As a first step in the proof of Proposition 13. 1[ note that since g has locally uni- 
formly controlled geometry and since v is uniformly bounded, we obtain directly 
from classical Holder estimates that v G A e (M). 



M&;0, 7 : = 
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Using this in (13.31) . we have 

((xd x ) 2 + A fc(0) - c{n)R m )v = -Ev + c(n)Av^ G A T \ (3.4) 

where r = 2±| (and we assume that re < 1 for simplicity). 

The conformal Laplacian Lg is an example of an elliptic 6-operator, and the 
operator on the left in this last equation is its asymptotic model at x = and called 
its indicial operator, I(Lg). This indicial operator can be analyzed quite directly 
using the Mellin transform in the x variable. To this end, introduce the indicial 
roots of L^; these are the values v for which there exists a function <p G C°°(Z) 
such that 

I{Lg)x v (t) = & L~ g x v (\) = 0{x v+t ') 

for some e' > 0. It is easy to see in this case, by separation of variables, that the 
indicial roots are given by 

v f := ± \Ai' where s P ec (-A(o)) = {Aj}- (3.5) 

The coefficient function (f> for any such indicial root equals the corresponding 
eigenfunction <pj. Since the lowest eigenvalue of — ^ = — A fc ( ) + c{n)Rk(o) 

is strictly positive, we have that ... < v± < Vq < < v q < < .... The 
indicial roots are the precise rates of growth or decay of approximate solutions of 

LgW = 0. 

The 6-calculus is merely a systematized method for passing from information 
about the indicial operator to the corresponding information about Lg itself. We 
quote some results from this theory, referring to ifTTl for a careful development of 
this 6-calculus as well as the more general edge calculus which will be invoked 
below. 

Proposition 3.3. ( IfTTl Theorem 4.4]) For k G N and < 7 < 1, the mapping 

Lg : x v C k b +2 >\M) — > x u C^'(M) 
is Fredholm if and only if v ^ vf for any j. 

Proposition 3.4. ( IfTTl Proposition 3.28]) Let f G A v ' (M) and suppose that 
LgV = /, where v G A U (M) for some v < v'. Then v G Aff g (M), or in 
other words, v has a partial expansion 

3=0 P =o 

where v G A v ' and the p,j lie in the interval (v, v'). Moreover, if f G Aff g (M), 
then v G A U A U '(M) and if f G then v G A" 
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Remark 3.5. Once we know that v G ^phg f or some v i v> > we can determine 
the exponents fij which appear in its expansion by a formal computation with the 
equation LgV = f. In particular, if the link metric k(x) depends smoothly on x, 
then all \ij are of the form Uj 1 + I, where is an indicial root and I G No- 

These are proved by constructing a parametrix G for Lg, which is a pseu- 
dodifferential operator, depending on the choice of (nonindicial!) weight v. The 
fundamental mapping results for this class of operators, proved in ifTTl §3] give that 

G:x^'- 2 ' 7 (M) -^x u Cl n {M), (3.6) 

G:A V \M) — >^(M), (3.7) 
G:A u phg (M) — ►.4£ hg (M) (3.8) 

are all bounded mappings, for any k G No, and v < v' with v ^ { v f}- 

Proof of Proposition \3. 1 \ Appealing directly to (13.41 ) and applying Proposition 13.41 
we deduce that v G -4phg- Using this on the right side of this equation gives 

2 

v G ^ph g e , an d bootstrapping further, we obtain that v G -4 p h g - The precise form 
of its expansion can then be determined by substituting this expansion into the 
equation. In particular, the leading exponent in the expansion is equal to one of the 
positive indicial roots , i.e. 

v~^(z)x u l +0{x u l +t '). 

Since v > when x > and the eigenfunction <pj changes sign unless j = 0, this 
expansion must start with 4>q(z)x p o . 

In the simplest case, where Rk(o) = (n— l)(n — 2), we can easily see that v§ = 
(n — 2)/2 and (f>o( z ) is a positive constant. In particular, u = x^ 2 ~ n ^ 2 v is bounded 
and strictly positive. In the more general case where we only assume that — ££> ) > 

0, we conclude that u ~ 4>o(z)x( 2 ~" n ^ 2+u o , where now 4>q is variable but still 
strictly positive. Recalling the conformal change g \-t g$ = x 2S ~ 2 g introduced 
in §2.3 which allows one to reduce to the case that u q = (n — 2)/2 with respect 
to the new radial coordinate £ = x s /S provided one chooses 5 correctly, we have 
that the solution is bounded and strictly positive relative to the background metric 
gs. Note that in the first case, the solution metric is still exact conic, while in the 
second case, it is conformally exact conic. □ 

3.2 Simple edges 

We next present the corresponding proof of regularity when (M,g) has simple 
edges. In this case we can take advantage of the construction of parametrices and 
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their mapping properties in the edge calculus. There is one important new step in 
the proof, beyond what was needed in the conic case, but then it proceeds exactly 
as before. 

We begin with exactly the same initial conformal change, replacing g by g = 
x~ 2 g, with the corresponding changes of conformal Laplacian and the Yamabe 
equation. We write the solution metric as v 4 ^ n ~ 2 ^g, so that the initial hypothesis 
is that < v < Cx € . Using the uniform geometry of g we deduce immediately 
that v is infinitely differentiable with respect to the geometry of this metric, or 
equivalently, in local coordinates, 

\(xd x y(xd y ) a dP v \ < c jja ,px e 

for all j, a, /?. Note however that this is not sufficient to assert conormality yet 
because these estimates do not control the y derivatives, i.e. the derivatives tangent 
to the singular stratum, as x — > 0. Obtaining this control requires the parametrix 
for the conformal Laplacian Lg as constructed in ifTTl . The difference between 
this operator and the one for the conic problem is the inclusion of the tangential 
part of the Laplacian x 2 A y ; thus when rewriting the Yamabe equation in the form 
(13 -4b . the error term E on the left includes x 2 A y v, which at this stage we cannot 
guarantee vanishes any faster than x e . 

The indicial roots of Lg are defined exactly as in the conic case above, and in 
particular are given by exactly the same formulas, with the important difference 
that the eigenvalues Xj of — can vary with y. Because of this, we do not 
expect solutions to have discrete asymptotics, i.e. polyhomogeneity, in the sense 
above, and we shall be satisfied with a more limited partial polyhomogeneity result. 
Recall, however, our key hypothesis that the lowest eigenvalue of — £^(o) * s e 1 ua l 
to c(n)/(/ — 1), and in particular is independent of y. Let i/ Q = ±^ c(n)f(f — 1) 
be the corresponding indicial roots. Fix < u < Vq. By lfl7l Theorem 6.1], there 
exists a parametrix G for Lg with the properties: G o L = Id — Q, where Q is 
a finite rank operator which maps into A"o +e ', and the analogues of (13.61 ), ( 13.71 ) 
and (13.81 ) are all valid provided v' < and furthermore, that we replace C^' 7 by 
the Holder spaces Ce ' 7 based on derivatives with respect to xd x , xd y and d z . The 
final extra property we need is that the commutator [d y ,G] = G\ is an operator in 
the edge calculus with exactly the same mapping properties as G itself; this is IfTTl 
Theorem 3.30]. 

We now use all of this information as follows. First let us apply the parametrix 

G to LgV + c(n)Az/ n+2 )/( n - 2 ) = 0, to obtain that 

n+2 

v = G(c(n)Av "- 2 ) + Qv . 
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We wish to use this equation to prove that v is conormal. Once this is done, the 
existence of the partial expansion follows directly from ( 13.7b and bootstrapping. 
Applying d y to each side gives 

n + 2 n+2 

d y v = G(c(n)kd y vn-->) + Gi(c(n)Av ™- 2 ) + d y Qv. 

The final term Qv is already conormal so presents no difficulties here. Apply- 
ing (13.61 ) to G\, we see that this term lies in x e Ce' 7 for all k. Finally, write 

Q yV (n+2)/(n-2) = cv V(n-2) d ^ and recall that w 4/(n-2) < Cx 4e/(n-2)_ j n ±Q 

simplest version of this argument, e > (n — 2)/4 so that v < Cx. Note that this 
is satisfied if we assume that our original solution u is bounded, which is a more 
natural assumption once we leave the setting of isolated singularities; for simplic- 
ity we assume that this is the case. Then |-y 4 /( n-2 ) d y v\ < \xd y v\, hence this term 
lies in r\^x e Ce ' 7 , and finally, d y v G C\kX e Ce ' 7 . Iterating this argument gives even- 
tually that <9° G D k x 6 Ce ' 7 for every multi-index a, so that u G A 6 . Proceeding as 
explained above shows that v, and hence u, has a partial expansion. 

3.3 The general case 

The final step in our proof of regularity is to extend these arguments to handle the 
case when M is a general smoothly stratified space and g is an iterated edge metric 
satisfying the hypotheses i) - iii) of §2.3. We also assume that the solution u is 
bounded. We shall be rather brief here since we have covered almost all of the 
main points of the argument already. 

Suppose that M is a space of depth k. By induction and the fact that this 
regularity theorem is localizable, we may assume that u has the appropriate regu- 
larity, i.e. partial polyhomogeneity, everywhere except possibly along the singular 
strata of highest depth, and so we can focus on these. Indeed, we can focus on the 
equation (13.11 ) in a neighbourhood of the form U x C\(Z) where U C M. y , Z is a 
compact smoothly stratified space of depth k — 1, where again the regularity result 
is known by induction, and C\(Z) is the truncated cone over Z. Let x be the radial 
variable on this cone. As before, we conformally transform this problem, writing 
the equation in terms of the new partially completed background metric g = x~ 2 g, 
so u = x( 2 ~ n ^ 2 v. As an important part of our inductive hypothesis, we assume 
the operator —C^^ on Z has discrete spectrum. This allows us to define indicial 
roots and analyze the indicial operator exactly as before. 

We use function spaces C*i' e 7 based on derivatives with respect to xd x , xd y and 
vector fields V which are tangent to the fibre Z and all of its singular strata. We 
have, by local elliptic regularity and induction, that v G x^ n ~ 2 ^ 2 C^ for all k, so 
the remaining job is to prove that v is conormal, that every tangential derivative 
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dy decays at this same rate. We do this by applying exactly the same commutator 
argument, but unfortunately there is no ready-made iterated edge calculus to which 
we can appeal. Fortunately we need very few consequences of such a calculus 
and can deduce these from a somewhat primitive parametrix construction. This is 
carried out in more detail in (3). Recall that we wish to construct an operator G 
such that, with L equal to the conformal Laplacian for g, GL = I — Q where Q 
maps into A U (M), which has mapping properties analogous to (I3.6I ). (I3.7I ). (I3.8I ). 
and finally, so that the commutator [d y , G] enjoys the same mapping properties. 

For simplicity let N denote the localized space U x C\[Z) blown up nbrewise 
at the vertex of each cone. Thus N is the product U x [0, 1) x Z. We construct G just 
as in [17 ] by regarding its Schwartz kernel G(x, y, z, x, y, 1) as a distribution on 
the space iV? e obtained from N x N by blowing up the fibre diagonal at the corner 
{x = x = 0}. In fact, this space is identified with the space (U x [0, 1))q X Z X Z, 
where the first factor is the 0-double space ofUx [0, 1), as constructed and used in 
irTTl . This has three boundary components, the left and right faces, corresponding 
to x — > and x — > 0, and the front face ff, which covers {x = x = 0, y = y} and 
is the face created in the blowup. The key point is that the lift of G to this space is 
conormal and partially polyhomogeneous at all faces. Its leading coefficient at ff 
is precisely the inverse for the so-called normal operator 

N(L) = (xd x ) 2 + x 2 A y + C n m 

which is globally defined on R + xl'x2. The invertibility of this normal operator 
on x v C^2 is tne main ellipticity hypothesis, and is proved exactly as in the simple 
edge case, relying on the fact that — ££( ) ^ s sa "i ct ly positive. We analyze this 
normal operator by taking the Fourier transform in y, thus reducing it to 

(xd x ) 2 + C n m -x 2 \r,\\ 

and then rescaling, setting t = x\q\, to arrive at 

(td t f + C n k{0) -t 2 , 

which is an operator on R + x Z. The inverse for this can be analyzed as in ifTTl . 
using mainly that ££(0) nas discrete spectrum. We denote the Schwartz kernel of 
this inverse by G(t, z, t, z). Rescaling and taking the inverse Fourier transform, we 
see that the Schwartz kernel of the inverse of the normal operator equals 

[ e^-y^Goixlvl^lril^,!)^ d7] 

for an appropriately chosen q. 
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The mapping properties for this inverse are deduced exactly as in the simple 
edge case. The final fact concerning the commutator [d y , G] can now be proved just 
as in the edge calculus, where this ultimately reduces to the fact that the commu- 
tator of the globally defined translation-invariant vector field d y on R £ commutes 
exactly with the explicit inverse for the normal operator written above. 

We have been (extremely) sketchy in the development in this last section. There 
are several reasons for this. The first is simply that while the idea is very close to 
that used in the simple edge case, it would still take considerable space to write out 
these details fully, and given the relatively minor importance of this final result, 
we have chosen not to do so. The sketch above is intended to provide a guide for 
anyone with a reasonable familiarity with the edge calculus. Finally, we point out 
that there are certainly other proofs that one might carry out to prove this regular- 
ity which would be more elementary in the sense that they do not explicitly use 
blowups and pseudodifferential operators, but which would require a substantial 
amount of verification of elementary details nonetheless. 

As described in the introduction, in a companion piece to this paper we give the 
full details of a proof of rather different sort of regularity statement which requires 
very little regularity of the background iterated edge metric, and which shows that 
the solution u to the Yamabe equation enjoys some Holder continuity properties. 

References 

[1] K. Akutagawa, Computations of the orbifold Yamabe invariant, Math. Z. 
271 (2012) 611-625. 

[2] K. Akutagawa and H. Kumura, Geometric relative Hardy inequalities and 
the discrete spectrum of Schrodinger operators on manifolds, to appear in 
Calc. Van 

[3] K. Akutagawa and B. Botvinnik, Yamabe metrics on cylindrical mani- 
folds, GAFA 13 (2003) 259-333. 

[4] P. Albin, E. Leichtnam, R. Mazzeo, P. Piazza, The signature package on 
Witt spaces, Annales de I'ENS 45 no. 2 (2012), 241-310. 

[5] B. Ammann, M. Dahl, E. Humbert, Smooth Yamabe invariant and 
surgery, arXiv:0804.1418, to appear in /. Diff. Geom. 

[6] B. Ammann, M. Dahl, E. Humbert, The conformal Yamabe constant of 
product manifolds, arXiv: 1103. 1826, to appear in Proc. Amer. Math. Soc. 



42 



[7] T. Aubin, Equations differentielles non lineaires et probleme de Yamabe 
concernant la courbure scalaire, /. Math. Pures Appl. 55 (1976) 269-296. 

[8] A. Bjorn and N. Marola, Moser iteration for (quasi)minimizers on metric 
spaces, Manuscripta Math. 121 (2006) 339-366. 

[9] H. Brezis and E.H. Lieb, A relation between pointwise convergence of 
functions and convergence of functionals, Proc. Amer. Math. Soc. 88 no. 
3, (1983) 486^90. 

[10] G. Carron, Inegalites de Hardy sur les varietes Riemanniennes non- 
compactes. /. Math. Pures Appl. 76 (1997) 883-891. 

[11] J. Cheeger, Spectral geometry of singular Riemannian spaces, J. Diff. 
Geom. 18 no. 4 (1983) 575-657. 

[12] T. Coulhon : Iteration de Moser et estimation Gaussienne du noyau de la 
chaleur, /. Operator Theory 29 (1993) 157-165. 

[13] D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of 
Second Order Grundlehr. Math. Wiss., vol. 224, Springer, Berlin, 1977. 

[14] M. Gursky, Compactness of conformal metrics with integral bounds on 
curvature, Duke Math. J. 72 no. 2 (1993) 339-3671. 

[15] T. Jeffres and J. Rowlett, Conformal deformations of conic metrics to 
constant scalar curvature, Math. Res. Let. 17 no. 3 (2010) 449^65. 

[16] J. Lee and T. Parker, The Yamabe problem, Bull. Amer. Math. Soc. 17 
(1987) 37-81. 

[17] R. Mazzeo, Elliptic theory of differential edge operators I, Comm. P.D.E. 
16 no. 10(1991) 1615-1664. 

[18] R. Mazzeo, Regularity of the singular Yamabe problem, Indiana Univ. 
Math. J. 40 (1991) 1277-1299. 

[19] R. Mazzeo, Resolution blowups, spectral convergence and quasi- 
asymptotically conical spaces. Journees Equations aux Derivees Par- 
tielles, Evian-les-Bains (2006) Expose VIII (available through num- 
dam.org). 

[20] R. Mazzeo and R. McOwen, Singular Sturm-Liouville theory on mani- 
folds, J. Diff. Eqns. 176 no. 2 (2001) 387-444. 



43 



[21] J. Nash, Continuity of solutions of parabolic and elliptic equations. Amer. 
J. Math., 80 (1958) 931-954. 



[22] M. Obata, The conjectures on conformal transformations of Riemannian 
manifolds, J.Diff. Geom. 6 (1972) 247-258. 

[23] P. Petersen, Riemannian Geometry, 2nd ed. Graduate Texts in Mathemat- 
ics 171, Springer- Verlag, Berlin, 2006. 

[24] M. Pflaum, Analytic and geometric study of stratified spaces, Lecture 
Notes in Mathematics 1768, Springer- Verlag, Berlin, 2001. 

[25] M. Reed and B. Simon, Methods of Modern Mathematical Physics II: 
Fourier analysis, self-adjointness. Academic Press, New York-London, 
1975. 

[26] L. Saloff-Coste, Aspects ofSobolev Type Inequalities, London Mathemat- 
ical Society Lecture Notes, 289, Cambridge University Press, Cambridge 
(2002). 

[27] R. Schoen, Conformal deformation of a Riemannian metric to constant 
scalar curvature, /. Diff. Geom. 20 no. 2 (1984) 479-495. 

[28] R. Schoen and S.T Yau, Conformally flat manifolds, Kleinian groups and 
scalar curvature. Invent. Math. 92 no. 1 (1988) 47-71. 

[29] N. Trudinger, Remarks concerning the conformal deformation of Rieman- 
nian structures on compact manifolds, Ann. Scuola Norm. Sup. Pisa 22 
(1968) 265-274 . 

[30] N. Varopoulos, Hardy Littlewood theory for semigroups, J. Funct. Anal. 
63 no. 2 (1985) 240-260. 

[31] A. Verona, Stratified mappings - structures and triangulability, Lecture 
Notes in Mathematics 1102, Springer- Verlag, Berlin, 1984. 

[32] J. Viaclovsky, Monopole metrics and the orbifold Yamabe problem, An- 
nates de L'Institut Fourier 60 no. 7 (2010) 2503-2543. 

[33] H. Yamabe, On a deformation of Riemannian structures on compact man- 
ifolds, Osaka Journal of Mathematics 12 (1960) 21-37. 



44 



